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Abstract 

Context. Ultracompact X-ray binaries (UCXBs) typically consist of a white dwarf donor and a neutron star or black hole accretor. 
The evolution of UCXBs and very low mass ratio binaries in general is poorly understood. In particular, the dynamical behavior of an 
accretion disk extending to a large radius (relative to the orbit) is unclear. 

Aims. We investigate the evolution of UCXBs in order to learn for which mass ratios and accretor types these systems can exist, and 
if they do, what are their orbital and neutron star spin periods, mass transfer rates and evolutionary timescales. 
Mefhods. We compute tracks of a binary containing a Roche-lobe overflowing helium white dwarf in which mass transfer is driven by 
gravitational wave emission. For different assumptions concerning accretion disk behavior we calculate for which system parameters 
dynamical instability, thermal-viscous disk instability or the propeller effect emerge. The significance of these processes during the 
evolution of an UCXB is considered. 

Results. At the onset of mass transfer, the survival of the UCXB is determined by how efficiently the accretor can eject matter in the 
case of a super-Eddington mass transfer rate. At later times, the evolution of systems strongly depends on the binary's capacity to 
return angular momentum from the disk to the orbit. We find that this feedback mechanism most likely remains effective even at very 
low mass ratio. In the case of steady mass transfer, the propeller effect can stop accretion onto recycled neutron stars completely at 
a sufficiently low mass transfer rate, based on energy considerations. However, mass transfer will likely be non-steady because disk 
instability allows for accretion of some of the transferred matter. Together, the propeller effect and disk instability cause the low mass 
ratio UCXBs to be visible a small fraction of the time at most, thereby explaining the lack of observations of such systems. 
Conclusions. Most likely UCXBs avoid late-time dynamically unstable mass loss from the donor and continue to evolve as the age 
of the Universe allows. This implies the existence of a large population of low mass ratio binaries with orbital periods ~ 70 - 80 min, 
unless some other mechanism has destroyed these binaries. Even though none have been discovered yet, black hole UCXBs could 
also exist, at orbital periods of typically 100 - 1 10 min. 
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1. Introduction 

Ultracompact X-ray binaries (UCXBs) are a subclass of low- 
mass X-ray binaries and consist of a white dwarf or helium star 
losing mass to a neutron star or black hole, at a sub-hour orbital 
period (S avonije et al.|1986 ). The short orbital periods point to 
white dwarf or helium burning star (sdB star) donors, as those 
are the only known stellar types that have the same size as the 
donor Roche lobe corresponding to these orbital periods ( Nelson 
|et al . igsF). The observed helium and carbon-oxygen composi- 
tion of the transfeiTed matter confirms this picture (Schulz et al. 

[2oor 



Nelemans et al. 2004 2006[ l. The accretor type can be 



identified via its inferred mass, magnetic field, spin period, ex- 
plosive nuclear fusion on its surface or inner accretion disk be- 
havior, and can be either a neutron star or a black hole (the latter 
have not yet been observed in UCXBs). To date about 30 UCXBs 
and candidates have been identified ( |in't Zand et al.|2007] l. 

UCXBs are important objects to study because of the ab- 
sence of hydrogen in the accretion disk and X-ray bursts (|in't| 
Zand et al. 2005[ l, their gravitational wave signal (Nelemans] 
2009|l, and them being te sts for the common-envelope phase 



( [Nelemans & Jonker|2010| ). Furthermore, UCXBs are candidate 
progenitors of radio millisecond pulsars ( [Alpar et al.||1982 1 be- 
cause of the high number of accreting millisecond pulsars found 



in UCXBs ( | Wijnands]|20 1 | . The ratio between the number of 
black hole and neutron star accretors constrains the high-mass 
end of the binary initial mass function. 

1.1. Evolutionary history 

The formation of ultracompact X-ray binaries starting from a 
zero-age main sequence binary involves one or two common en- 
velope stages and a supernova. The initially more massive star 
evolves off the main sequence and starts losing mass to its com- 
panion. Depending on the mass ratio and stage of donor evolu- 
tion, this may happen in an unstable way, leading to a common 
envelope. The orbit decays significantly due to energy and angu- 
lar momentum loss via friction. If the system does not merge, it 
emerges as a much shorter-period system consisting of the core 
of the giant and a rejuvenated main sequence companion. There 
are several scenarios for the subsequent evolution; the core of 
the giant may develop into a core-collapse supernova, or the 
main-sequence companion may evolve into a (sub)giant, caus- 
ing a second common envelope before it becomes a supernova 
itself ( [Tutukov & Y ungelso n|[T993| [Podsiadlowski et al.|[2002l l. 
If precisely one of the stars becomes a supernova and the sys- 
tem is not unbound during this event, the system can develop 
into a neutron star/black hole - white dwarf/helium star binary. If 
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neither of the components is massive enough to become a core- 
collapse supernova, a second common envelope can lead to a 
double white dwarf or a white dwarf - helium star binary. An 
accretion-induced collapse of a white dwarf may still yield an 
UCXB, once mass transfer has resumed ( |van den Heuvel|1984| l. 

More models for UCXB formation exist. In globular clus- 
ters, home to about one-third of the presently known UCXB 
population, formation may be dominated by dynamical inter- 
actions ( |Verbunt|[T987) [Ivanova et aL]|2005j |Voss & Gilfaiiovl 
2007| l. After a system with a neutron star or black hole compo- 



nent has formed, UCXBs can evolve through 3 scenarios: via a 
white dwarf donor (Yungel son et al.||2002[ |, helium-star donor 
(Savonije et al. 1986; , Yungelson 2008|l or a main-seque nce 



donor ([Podsiadlowski et al.|2002[|van der Sluys et al.|2005|pa 
|& Li|2009| l. However, after a relatively short time of mass trans- 
fer, each donor becomes degenerate and the subsequent evolu- 
tion is very similar for each scenario. Therefore, we consider 
only the white dwarf donor scenario. 

As suggested by the rarity of UCXBs, the initial stellar and 
binary parameters must be finely tuned in order to arrive at a 
stable ultracompact configuration. Several things can go wrong: 
binaries can either evolve too slowly to have two post-main se- 
quence components, be too wide to have interacting compo- 
nents, merge during a common envelope, possibly explode as 
a single degenerate type la supernova, become unbound by a su- 
pernova, or experience two supemovae. 

1.2. Present research 

The objective of this paper is to investigate the occurrence of 
a dynamical instability in low mass ratio UCXBs, and the in- 
fluence on UCXB evolution of magnetosphere-disk interactions 
and the thermal-viscous disk instability. By improving our un- 
derstanding of all evolutionary stages of UCXBs, we will be able 
to better estimate their observational properties, such as X-ray 
luminosity and orbital period, fraction of time visible and life- 
time. In a forthcoming paper, we will use the results to predict 
the present-day UCXB population, which can be compared to 
observations by the Galactic Bulge Survey ( Jonker et al.|201 1 1. 

In particular the behavior of old, low mass ratio (q < 0.01) 
and low mass transfer rate systems is unclear. None have been 
discovered, even though their long evolutionary timescales sug- 
gest that many should exist if no disruptive process emerges at 
some point. Are they invisible to our instruments (most of the 
time), or do they no longer exist? 

Section|2]gives an overview of regular UCXB evolution and 
the main angular momentum flows in the absence of instabilities. 
Several complications that may play an important role during the 
evolution of an UCXB are mentioned in Sect. |2.6| In the method. 
Sect. [3] we describe the evolutionary tracks as a function of the 
degree of feedback of angular momentum from accretion disk to 
orbit, the propeller effect and the disk instability model. In the 
results. Sect. [4] we analyze the tracks, consider how likely lim- 
ited feedback is, and look into the magnetic field, both for con- 
stant mass transfer and non-steady behavior such as outbursts. 
In Sect. [5] we discuss the results and give a conclusion. The ap- 
pendix contains fits of evolutionary tracks and some analytical 
approximations. 



2. Overview of UCXB orbital mechanics 

The evolution of a binary consisting of only degenerate stars or 
black holes cannot be driven by single star evolution, since such 



components hardly evolve, if at all, and certainly do not increase 
their radius in the absence of interaction with other stars. Since 
we assume white dwarfs to be cold, i.e. thermally relaxed, angu- 
lar momentum loss and redistribution are the only processes that 
can alter the system. 

2.1. Angular momentum flows 

The angular momentum flows in an UCXB are given by 



■^oib — -^1 



GWR 



■^stream ~^ >^torque ~^ >^eject5 



(1) 



where /orb < is the change in orbital angular momentum, 
■^GWR < the loss to gravitational wave emission, 7stream < an- 
gular momentum advected from the donor to the accretion disk 
along with the transferred matter, /torque > is the return of an- 
gular momentum from the disk to the orbit by means of a tidal 



torque between the oute r disk and the donor (Lin & Papaloizou 



1979 



Frank et al. 



2002 1 and /eject < is the angular momentum 
carried by matter ejected from the system. The high surface grav- 
ity of degenerate stars prohibits significant wind mass loss from 
the donor, so spin angular momentum loss via magnetic braking 
can be neglected. Change in accretor spin can present a sink or 
source of angular momentum ( [Marsh et al.||2004| ), but contrary 
to white dwarf accretors, this is negligible compared to the an- 
gular momentum transport to and from the disk (Priedhorsky & 
|Verbunt|l988) for neutron star and black hole accretors. 

When the system has an accretion disk, 7stieam is precisely 
balanced by the feedback flow in the opposite direction, /torque 
(Priedhorsky & Verbunt 1988). In that case, the accretion disk 
does not gain or lose angular momentum, which leaves orbital 
angular momentum loss via gravitational wave radiation as the 
sole process driving binary evolution, giving the usual expres- 
sion 



J, 



orb 



J, 



GWR 



+ J, 



eject- 



(2) 



However, below we will consider the more general case. The 
orbital period of two stars must be rather short for gravitational 
wave emission to be able to significantly shrink the orbit within 
the age of the Universe, e.g. shorter than 14 hr in the case of two 
1 Mq components. Once the orbit has shrunk to the point where 
the Roche lobe of a white dwarf component is smaller than the 
associated star's volume, the white dwarf starts to lose surface 
matter 

2.2. Stable mass transfer 

To calculate the mass transfer rate, we assume that the donor 
precisely fills its Roche lobe at all times t during stable mass 
loss (in reality the donor slightly overfills its Roche lobe, but this 
distance is small compared to the Roche-lobe radius and can be 
neglected): 



Rd(t)^RL(t) and Ra(t) ^ R^t) Vf 



(3) 



where R^ is the donor radius and 7?l is the effective radius of the 
donor Roche lobe. Equation |3]implies that the ratio of the com- 
ponent masses and the donor radius determines the semi-major 
axis and orbital period, which increase with time. The change 
in R^ and R]^ can be separated in a mass-transfer dependent and 
-independent part as ( [Soberman et al.|I997| l 
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Figure 1. Tracks for a helium white dwarf donor with an initially 
1.4 Mq, 12 km radius neutron star (solid) and 10 Mq black hole 
(dashed) accretor The tracks start at the respective dynamical 
instability limits of 0.83 Mq and 1.04 Mq, since systems can 
only survive if they start with a lower donor mass (Sect. 2.3 1. 
The isotropic re-emission limit is ignored in the first part of the 
track. The numbers correspond to the two circles on their right- 
hand side and represent the logarithm of the system age (yr), 
which is the time since the onset of mass transfer. The relative 
location of both tracks will be discussed in Sect. 14.2.21 



as it loses mass, and since it is assumed to fill its Roche lobe 
at all times (Eq. |3]l, the orbital separation and hence orbital pe- 
riod must increase as well. Since gravitational wave radiation 
becomes significantly weaker in a wider orbit, the mass transfer 
rate decreases. This can be understood as follows: mass transfer 
from a less massive to a more massive component in itself tends 
to widen the orbit, and the weaker the orbit-shrinking gravita- 
tional waves, the less mass transfer is needed to keep the orbit 
wide enough to accommodate the donor 

If the binary transfers too much matter, the donor detaches 
and mass transfer stops. Conversely, if too little matter is trans- 
ferred, the donor overfills its Roche lobe more and more, increas- 
ing mass loss. An equilibrium is reached; the system naturally 
arrives at a mass transfer rate for which the donor overfills its 
Roche lobe by the right amount. 

2.3. Dynamical stability 

From Eq. (|6| it can be seen that the mass transfer rate becomes 
unbounded when and approach each other Dynamical sta- 
bility of the donor requires that^pon mass loss, the change in 
/?L exceeds the change in R^, i.e|^ 



Ri. > Rd. 



(7) 



Dividing both sides by R^ - R j^ (Eq.[3} a nd by Md/M^ < yields 
the well known criterion (e.g. Webbink ( 1985[ l; Pols & Marinus 
(fl994|f] 



(8) 
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where (\ = d\nRJd\nMi (Hut & Paczynski 1984 1 represents 
the change in R{ resulting from mass transfer or loss, i.e. a 
change in donor mass (Ri oc M^'). The subscript evo denotes 
single star evolution and ami denotes angular momentum loss 
from the system (in the absence of mass transfer and loss). The 
advantage of the above separation of the change in radii is that 
it allows for more easily solving the equilibrium mass loss rate. 
From Roche lobe geometry and the orbital angular momentum 
equation follows that (d In ^L)ami = (d In fl)ami = 2(5 In yorb)aini in 
the case of fixed component masses (a is the semi-major axis). 
Combining the above and taking (51n/?d)evo - for the fully 
cooled, non-evolving, white dwarf case, leads to the (negative) 
donor mass loss ratqj 



Md 
Md 



id-(i 



(6) 



where orbital angular momentum (unrelated to mass transfer or 
mass loss) is lost exclusively by gravitational wave emission (Eq. 
|2]), given by ( Landau & Lifshitz||1975 i. The eff'ect of angular 
momentum loss via mass ejection is contained in and will be 
discussed in Sect. |3.2| In Fig. [T] we show the mass transfer rate 
given by Eq. (j6|l for ^ given in Sect. 3.3 



Figure[T]shows how the mass transfer rate decreases with in- 
creasing orbital period and time. The white dwarf grows in size 

' An approximate analytic solution of Eq. ^ is given in Appendix 
[B| along with fitted tracks in Appendix [A| 



Physically this means that when a donor that does not obey this 
criterion starts transferring matter via Roche-lobe overflow, its 
volume grows faster than the volume of the Roche lobe it is con- 
tained in. The donor will overflow its Roche lobe more and more, 
leading to runaway mass loss on the dynamical timescale of the 
donor (~ 5 s for an 0.6 Mq white dwarf), and followed by dis- 
ruption of the donor 

Near the dynamical instability limit, the mass transfer rate 
becomes in principle arbitrarily high, as shown by the upper left 
part of the tracks in Fig. [T] This means that in order to find the 
critical mass ratio at which this instability occurs, we have to 
consider for the case in which a fraction 1 of transferred mass 
is ejected. Figure[2]shows that for a system with an accretor mass 
Ma = 1 .4 Mq (dotted curve), mass transfer is dynamically unsta- 
ble if Md > 0.83 Mq. In the case of a 10 Mq accretor (dash- 
dotted) this becomes Md > 1.04 Mq. 

2.4. Isotropic re-emission 

Even when the donor mass is low enough to avoid a dynamical 
instability, the high amount of mass it needs to lose to avoid a 
merger of the binary components shortly after the onset of mass 
transfer ( [van den Heuvel & Bonsema, 1984) may be too much 
for the accretor to accrete or eject. The ultimate condition for 



- Note that this criterion has a different meaning than the second part 
of Eq. ([3]l. Here, change in radii resulting from only mass loss is con- 
sidered, whereas Eq. ^ applies when gravitational wave radiation is 
included as well. 

^ A stability criterion expressed in ^-values, i.e. logarithmic deriva- 
tives, is to be preferred over linear derivatives of radius to donor mass 
because the expressions for and ^j, Eqs. l |13| ) and \26\ respectively, 
are simpler than the expression for dRJdMi = <fi • RJM^ as the factors 
Ri and introduce several extra variables. 
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Figure 2. Mass-radius relation exponent (Sect. 3.2 1 for a 
1.4 Mq accretor (dotted) and a 10 Mq accretor (dash-dotted) and 



the zero-temperature white dwarf (dashed. Sect. 3.3 1. Nearly 
all matter is assumed to be lost via the isotropic re-emission 
mechanism (Sect.[Z4]i. 



survival is how efficiently the accretor can eject matter in the 
case of a super-Eddington mass transfer rate. The upper limit 
of mass that can be removed from the system follows from the 
energy balance. By assuming that precisely the Eddington limit 
is accreted, the large amount of gravitational energy liberated 
can in principle be employed to unbind additional transferred 
matter (with energy corresponding to the first Lagrangian point) 
from the system. This mechanism is called isotropic re-emission. 
Soberman et al.|(| 19971 and Tauris & Savonije (1999) described 
the mechanism as matter transferred via an accretion disk to the 
vicinity of the accretor, from where it is subsequently ejected as 
a fast, isotropic wind. 



2.5. The deciding stability criterion at the onset of mass 
transfer 

The onset of Roche-lobe overflow is a very important event in 
the evolution of a binary. White dwarf donors that are too mas- 
sive lose mass in a dynamically unstable way (Sect. [23] l and the 
systems containing these do not survive. Figure [3] (upper solid 
curve) shows that systems with a 1.4 Mq neutron star accretor 
are dynamically unstable if the white dwarf donor is more mas- 
sive than 0.83 Mq. Since there is no accretion disk present yet 
when mass transfer first starts, there is no tidal torque transfer- 
ring back angular momentum to the orbit (to be discussed in 
Sect. 2.6. l| l, and the criterion becomes stricter; the system must 
lie within the lower solid curve in Fig. [3] However, if the donor 
is less massive than ~ 0.50 Mq (dotted curve), the disk forms on 
a timescale short enough to start transferring back enough an- 
gular momentum to stop the dynamical instability (|Verbunt^&j 
Rappaport||1988 1. Moreover, the 0.50 Mq disk formation crite- 
'aced by an even stricter criterion: the isotropic 
0.37 Mq (dashed curve)01n this pa- 



rion itself is rep 
re-emission limit of Mj 

per we assume that isotropic re-emission is the criterion which in 



[Yungelson et al.| ( |2002| l found a value of 0.44 Mq, for a 1.433 Mq 
neutron star accretor, due to a different adopted isotropic re-emission 
prescription combined with the relatively weak sensitivity to donor 
mass of both mass transfer rate and isotropic re-emission rate. 



Dynamically stable with feedback 




4 6 
Accretor mass (Mq) 

Figures. Stability regimes in terms of zero-temperature helium 
white dwarf donor mass versus accretor mass. The upper solid 
curve is the upper limit for dynamical stability in the presence of 
an accretion disk that feeds back angular momentum, where all 
transferred matter is lost via isotropic re-emission. The dashed 
curve is the upper limit for isotropic re-emission (i.e. if the 
system lies above this curve, the accretor will be engulfed by 
transferred matter). The lower solid curve encloses the parame- 
ter space where systems are dynamically stable in the absence 
of angular momentum feedback. The dotted curve indicates a 
very rough estimate of the upper donor mass for which there 
is enough time to form an accretion disk and suppress the no- 
feedback instability, based on jVerbunt & Rappaport ^1988) . 



practice determines whether a system survives the onset of mass 
transfer, or merges]^ 

If the binary system survives the onset of mass transfer, it 
will settle in a state of stable mass transfer, as described in Sect. 



2.6. Complications 

The perpetually expanding binary with declining mass transfer 
rate described above is idealized. In reality, several processes can 
potentially disturb continuous accretion or even disrupt the bi- 
nary, and they may become important in view of the low mass ra- 
tios and mass transfer rates expected for compact systems driven 
by gravitational wave radiation. 

2.6.1 . Limited feedback of angular momentum at low mass 
ratio 



In Sect. 2.3 the possibility of a dynamical instability at the onset 
of mass transfer was discussed. According to Fig.|2]such an in- 
stability can only occur at a high donor mass of M^j > 0.8 Mq, 
depending on accretor mass. This is no longer necessarily true 
if an extra angular momentum sink is introduced. Here we look 
into possible reasons for limited feedback of angular momentum 
from disk to orbit. 



^ If isotropic re-emission would fail to prevent a merger, UCXB could 
only be formed via an initially non-degenerate donor like a helium star, 
with the constraint that it should not become fully degenerate before 
reaching a mass below 0.08 Mq, when mass transfer becomes sub- 
Eddington. 
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Angular momentum is advected from the donor to the ac- 
cretion disk along with transferred matter. Initially, when a disk 
has yet to form, matter accumulates in a ring around the accretor 
and, due to viscous friction, spreads out radially to form an ac- 
cretion disk. The disk keeps expanding until it reaches a radius 
that is so large that all added angular momentum can be returned 
to the orbit via the tidal torque between the donor and the outer 
disk, and no further (see | Frank et al.||2002 | chap. 5). A torque 
exists because the disk is slightly elongated due to the tidal force 
of the donor, and also slightly ahead of the donor because of the 
shorter orbital period and a delay in adjusting its shape. At the 
same time, the outwards transport of angular momentum through 
the disk allows for matter to move inwards, and eventually to be 
accreted onto the central object. 

Possible causes for reduced feedback include a viscosity that 



is too low and tidal torque that is too weak (Lin & Papaloizou 
[1979 ; Ruder man & Shaham 1983, H ut & Paczynski|1984| l, and 
the resumption of mass transfer after a detached stage caused 
by orbital widening, which in turn is due to mass loss caused 
by an accretion-induced collapse of a white dwarf accretor ( van 
den Heuvel 1984). Recently, [Yungelson et al.| ( |2006| l; 



Lasota 



( |2007| l stated that it remains unknown what happens to mass 
transfer and the accretion disk in systems with a mass ratio 
q — Md/M^ < 0.01. When q < 0.02, the circularization radius 
exceeds the estimates of the outer radius by Paczyhski ( 1977| l 
and Papaloizou & Pringle ( 1977} , causing the disk to be trun- 
cated, and matter to circularize onto unstable orbits, thereby 
preventing accretion (l Yungelson et al.|2006[|Lasota|2008| l. The 
donor may cause gaps in the disk at certain radii via destabiliz- 
ing orbital resonances, causing disk particles to be confined to 
certain radius ranges. This process is similar to the formation of 
gaps in the ring system of Saturn by Mimas and other moons 
(Franklin & Colombo 19701), and the Kirkwood gaps in the as- 



teroid belt due to Jupiter's gravitational influence. Matter may 
be unable to pass these gaps, which would prevent both the out- 
ward expansion of the disk to the radius required for full angular 
momentum feedback to the orbit, as well as the corresponding 
accretion on the inside. 



2.6.2. Occurrence of a dynamical instability at low mass ratio 

A dynamical instability at low mass ratio may happen only if 
feedb ack would be reduced for any of the reasons given in Sect. 



2.6.1 



In the absence of feedback ytorque = and Eq. (|Tj) becomes 



Joih - JgWR + Jsl 



ejects 



(9) 



The disk would turn into an extra angular momentum sink and 
the orbit loses more angular momentum than in the case of feed- 
back, and has a stronger tendency to shrink. More mass must be 
transferred to compensate for this and keep the donor within its 
Roche lobe, and thereby a new equilibrium is established. For 
systems with a low mass ratio this effect is stronger because of 
the geometry of the system: mass entering the disk carries a high 
specific angular momentum because of the large distance (rela- 
tive to the semi-major axis) of the first Lagrangian point to the 
center of mass. At very low mass ratio, the orbit loses so much 
angular momentum to the disk that no amount of mass transfer 
can compensate for this; the donor will be disrupted]^ 



* The dynamical instability for high mass white dwarf donors on the 
other hand is a result of the steep mass-radius relation of the donor, and 
also because the orbit expands less easily when the mass ratio is closer 
to 1 (Fig.|2j. 



[Ruderman & Shaham| ( T983]) wer e the first to identify this 
instability, and Ruderman & Shaham ( 1985") showed that in the 
case of a helium white dwarf donor, the dynamical instability can 
be reached after 9 Gyr, with Mj-A - 10"^ Mq for M^- 1.4 Mq. 



Bonsema & van den Heuvel ( 1985 1 and [Verbunt & Rappaport 
( |1988| l came to a similar conclusion . 



On the other hand, "Jeffrey ( 1986| l reported that the fast spin 
period of the single millisecond pulsar PSR 1937+214 contra- 
dicts a history in a binary followed by donor disruption, because 
the spin would have been much slower at the evolutionary stage 
such a disruption would take place. 

Bildsten (2002) considered finite-temperature white dwarfs 
and found that those would reach the instability at higher donor 
mass than zero-temperature white dwarfs, for instance at Mj = 
0.01 in the case of a 10^ K helium white dwarf donor and a 
Ma = 1 .4 Mq accretor. 

The questions we want to answer are, 1) could feedback be 
reduced at some stage in the evolution of an UCXB, and 2) if 
so, does this lead to a dynamical instability within the age of the 
Universe? 

2.6.3. Propeller effect 

More processes can potentially limit accretion. In the case of a 
neutron star accretor, accretion may be disturbed by the mag- 
netic field. This is the propeller effect, in which the fast-rotating 
magnetosphere accelerates transferred matter in the azimuthal 



direction and stops it from accreting (Davidson & Ostriker 1973 
Illarionov & Sunyaev,1975j . 



2.6.4. Thermal-viscous disk instability 

At low mass transfer rates, the disk instability model describes a 
thermal-viscous instability resulting from a relatively large and 
sudden local increase in opacity and viscosity of the disk ma- 
terial (Osaki 1974 Lasota 200 Ij ) (note, this disk instability is 
unrelated to the dynamical instability discussed in Sects. |2.3| 
and 2.6.2 1. The high viscosity causes a much higher accretion 
rate, called an outburst. Outbursts are alternated by low-viscosity 
stages during which the disk builds up again. Stable behavior 
persists if the entire disk has a rather homogeneous degree of 
ionization; UCXBs in particular can have a stable disk if the 
mass transfer rate is sufficiently high to keep the entire disk ion- 
ized through X-ray irradiation ( in't Zand et al.|20O7] l. 



3. Method 

In this section we describe the details of our treatment of the 
evolutionary stages of UCXBs. The concepts that will play a 
role during the evolution of UCXBs are discussed in order of 
appearance. 

3.1. Isotropic re-emission 



As discussed in Sect. 2.4 the energy released by accreting some 
matter, can unbind additional matter via radiation pressure. The 
highest spherical accretion rate possible without the arriving 
matter stopped and blown outwards by radiation pressure is the 
Eddington accretion limit Mndd = ^^cR^/Kss, with c the speed of 
light, /?a the accretor radius and /Ces ~ 0.199 cm^g"' the elec- 
tron Thomson scattering mean opacity for hydrogen-deficient 
matter The maximum mass transfer rate that can be survived 
via isotropic re-emission is given by the ratio of the potential 
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well at the accretor surface Oa and the potential in the first 



Lagrangian point (Begelman 1979 King & Begelman 1999 



Yungelson et al. 2002 1. The potential in LI by definition is equa 
to the corotating-binary potential at the Roche lobe surface Ol, 
which can be approximated to within 0.15% by 



GMtot/a 2 



= T + exp < - 



[in(lO) 



12 



9/&\ 



7/16) 



(10) 



with < ^ < oo, Mtot = Ma + Md the total binary mass and G the 
gravitational constant. If we define (piq) as the right-hand side of 



Eq. ( 10 1, the maximum survivable mass transfer rate is 

■GMJR., 



4nc 



-(f)(q)GMtot/a A-es (/)(q)(l+q) 



(11) 



and is independent of R^. Since the matter is ejected from the 
vicini ty of the accretor ( [Soberman et al.|1997{|Tauris & Savonije] 
1 1999 1), it has approximately the same specific orbital angular mo- 
mentum as the accretor. 



3.2. Response of Roche-lobe radius to mass transfer 

Any effect of change in angular momentum because of mass 
transfer and ejection from the system is contained in - 
(91n/?L/51nMd (Sect. |2.2|i. The eff'ective radius of the donor 



Roche lobe /?l can be approximated to within 1 % by Eggleton 
( [T983l l 



0.49 



0.6-H^-2/3ln(l 



(12) 



Given q = M^/M^ and R^, Eqs. Q (first part) and ( [T2] l yield a. 
The Roche-lobe geometry only applies when the spin period of 
the Roche-lobe filling component is equal to the orbital period, 
and when the orbit is circular. These conditions are normally met 
in compact Roche-lobe filling systems. 
Logarithmic diff'erentiation of yields 



where = dlna/dlnM^ and 



(13) 



with 



dln{Ri^/a) 
dlnMd 



21n(l +^i/^)-(l +q-^'^y 



ln(l -H ^1/3) H- 0.6^2/3 



= d\nq/d\nMa = 1 -H (1 - /3)q, 



(14) 



(15) 



where /? is the fraction of mass lost by the donor that leaves the 

^ with /J 



Ois 



system via isotropic re-emission, see Eq. (21 
the same as in Ma rsh et aL] ( |2004l ), who stuHied mass transfer in 
binary white dwarfs. 

3.2. L Mass flows 

In order to determine we consider three mass flows in the 
binary system: mass lost by the donor can either be accreted, 
stay in the disk or be unbound from the system. 




0.4 0.6 0.8 

White dwarf donor mass (Mq) 



Figure 4. Accretion efficiency e (Eq. 17 i for a white dwarf donor 
and an initially 1.4 Mq neutron star (solid) and a 10 Mq black 
hole (dashed) accretor, in the case of full feedback of angular 
momentum. Both curves start at their respective dynamical sta 
bihty Hmit (Sect 



2.3 1 



Accreted mass The accretion efliciency e (where < e < 1) is 
the fraction of the matter lost by the donor that is accreted by the 
primary 



Ma 



-eMd. 



(16) 



When we assume full accretion in the case of a sub-Eddington 
mass transfer rate (in the absence of magnetic fields or instabili- 
ties) 



. , MEdd 

e = nun ( — , 1). 

Md 



(17) 



Medd ~ 3 ■ 10" MQyr" for a neutron star that is accreting 
hydrogen-deficient matter The evolution of e for two different 
accretor masses is illustrated in Fig.|4] 

Mass lost from the system In the case of non-conservative 
mass transfer the system loses angular momentum along with 
the ejected matter The specific angular momenta of the ejected 
matter and the binary are 7oib/Mtot and 7oib/Mtot, respectively. 
The parameter y is defined as the ratio between these 



(j/n 



orb 



M,o,/M,c 



(18) 



(note that y in |Soberman et al.| ( |1997| has a different meaning, 
unlike and e). It can easily be derived that y = ^ in the case 
of isotropic re-emission, because matter has the specific orbital 
angular momentum of the accretor (Sect.jTTJ. 

Mass added to accretion disk In the absence of an accre- 
tion disk that can efficiently transfer angular momentum back 
to the orbit, the stream and torque terms in Eq. ([T|l do not can- 
cel out, and the disk becomes a net sink of angular momentum. 
In the extreme case of no feedback at all /torque - 0, leaving 



-/stream = VCM^^Md ( jllut & Paczyhski 1984]!. The circulariza- 
tion radius Rh is the orbital distance at which the specific angular 



6 



van Haaften et al.: The evolution of ultracompact X-ray binaries 



momentum is the same as in the incoming stream from the donor 
(Osaki 1989 1, minus the amount of angular momentum removed 
by the donor, approximated by ( ^Verbunt & Rappaport,1988j 

= 0.0883 - 0.04858 logio(^) 

+0.11489 logJ(,(^) + 0.020475 log^o(^), 

adjusted for our definition of ^ = M^/M^, and valid for 10"^ < 
q < I. 

In the general case of limited feedback, we introduce a factor 
J] (where < ?; < 1) describing the fraction of mass lost from 
the donor that is (net) added to the disk because it cannot be 
accreted, defined analogous to e in Eq. ([T6]l 



(19) 



(20) 



Conservation of mass is represented hy j3 + e + t] = l.ln the 
case of full feedback (j] = 0), /3 may be estimated by assuming 
that all non-accreted matter leaves the system, unless the maxi- 



mum isotropic re-emission rate M„ 
accomplish this, so 

y6 = mm( — 1 - e). 



> (Eq. 



11 



is too low to 



(21) 



3.2.2. Response of semi-major axis to mass transfer 

Combining the above equations, the change in angular momen- 
tum of the system j^^-Q because of mass leaving the system or 
entering the disk is 



J ^ Tl 

Miol J orb 

Mtot rr, 



(22) 



with - Rh/a. 

To find (a, we first differentiate the logarithm of the equation 
for orbital angular momentum to time. We introduce a disk mass 
term to allow for transferred matter to stay in the disk instead 
of being accreted, corresponding to possible limited feedback to 
find 



orb 



_ 1 fl Md Ma + Mdisk 1 Mtc 

~ 2 fl Md Ma 2 M,c 



(23) 



where Mdisk ^ Mg has been used. Mass added to a disk around 
the accretor has the same effect on Jorh as mass that is accreted, 
because the accretor-disk subsystem is orbiting the center of 
mass of the binary. Hence the generalization of accretor mass 
to include disk mass. 

Equating Eqs. (22 1 and (23 i eliminates (7/7)oi-b. Inserting 



Eq. ( 16 1, Mtot = ySMd and M^ + Mdisk = (P - l)Md gives 



l-(l-/3)q- 



1+^-1 



(24) 



In the case of full feedback, 77 = and [i + e - 1 . At the onset 
of mass transfer from a massive donor, when a disk has not yet 
formed, rj - \ and - -2[1 - q - -^(1 + q)rb\, identical to the 
equation in |Verbunt & Rappaport| ( |1988| l; |Marsh et al.|p004| l. In 



the case of low-mass donors where accretion is made impossible 
by limited feedback, accretion has ceased (e - 0) and isotropic 
re-emission has stopped accordingly (/? = 0), hence rj - \. The 
actual value of 77 during the evolution of an UCXB will be con- 
sidered in Sect. 14. II 

The relation as a function of q for e = is shown in Fig.|2] 
for two accretor masses and for more general mass flow param- 
eters in Fig.|7] for only a neutron star accretor The Roche-lobe 
radius depends on both component masses and angular momen- 
tum losses via its dependence on a, and on the mass ratio via 

Note that the variables e, Md, ^l, (-//-/)gwr, Md, Ma and 
Mfidd, and in the case of a black hole accretor also R.^, are cir- 
cularly related to each other through Eq. (j6]l, the gravitational 
wave equation, and Eqs. ( 17 1 and ( [T3] l. Hence, iteration at each 
mass step is necessary to arrive at the equilibrium values. 

3.3. White dwarf radius 

Now is known, only the response of the white dwarf donor 
radius to mass loss, ^d, is still needed to evaluate Eq. (j6]). We 
justify the use of a zero-temperature white dwarf mass-radius 
relation by noting that thermal pressure is negligible for sig- 
nificantly degenerate objects such as white dwarfs with a mass 
above 0.01 M©, and still quite small for lower masses. From the 
mass-radius relations derived by Deloye & Bildsten^ (2003) it 
follows that the effect of ideal gas pressure on the radius, ad- 
ditional to degeneracy pressure and Coulomb attraction, is less 
than 10% for Md > 10"^ Mq when the temperature T - l(f K, 
less than 2.3% for Md > 4 ■ 10"^ Mq when T = 10^ K and less 
than 0.1% for Md > 10"^ Mq when T = 10"^ K. Temperature 
is important only for very low-mass or very hot donors. For a 
mass ratio below ~ 0.01, donor radii could be 10 - 30% higher 
than for the zero-temperature case ( [Nelso n & Rappaport|2003j. 
Furthermore, we consider only degenerate donors, since any fu- 
sion will be extinguished on a relatively short timescale after 
some mass has been lost ( Savonije et al.|1986 Yungelson 2008| l. 
Also, a non-degenerate outer layer can have a significant effect 
on the radius response to mass loss, but again this is temporary 
- once the layer has been lost, a degenerate object remains. 

The radius R^ of a zero-temperature white dwarf (the ratio 
of atomic number to atomic weight Z/A - 1/2) of mass Md by 
Eggleton ( Rappaport et al.|19^7] l is 



Ro 



0.0114 



\Maj 



-2/3 



2/3 



1/2 



2\mJ \mJ 



-2/3 



(25) 



with Moi =1.44 Mq. For a helium white dwarf, Mp - 5.66 ■ 
10"* Mq and for a carbon-oxygen white dwarf (assuming 50% 
of ea ch element) Mp - 1.44 ■ 10"^^ Mq, valid for < Md < 
Mch ( |Marsh et al.|2004| l. ^d = dlnRa/dlnMa is determined by 
logarithmic differentiation 



1 + 



4/3 



( 



Ma) 



1 + 



I \l/3 
3 l^Mpj 



1 + 



(26) 



and is illustrated by the dashed curve in Fig.|2] Throughout this 
paper we consider helium white dwarfs, though all results except 
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quantitative details and the disk instability details also apply to 
carbon-oxygen composition. 

3.4. Propeller effect 

3.4.1 . Magnetic field of neutron star accretors 

Neutron stars typically possess a strong magnetic field. During 
accretion, the magnetic field strength is reduced ( Bhattacharya| 
& van den Heuvel„1991j , but even for old millisecond pulsars 
with a spin period below 30 ms, it can still be 10^"^ G (Wang 



|et aLII^Ol T). The magnetic field which is locked to the neutron 
star is supposed to penetrate the accretion disk, and in the case of 
a rapidly spinning neutron star, the outer magnetosphere forces 
the orbiting matter to velocities exceeding the local Kepler ve- 
locity by exerting a torque. As a result of the angular momentum 
gained, disk material in this region moves outwards, which re- 
duces accretion and can lead to unbinding the matter from the 
system. 



3.4.2. Critical mass transfer rate 

The Alfven radius (also known as magnetosphere radius) is the 
distance from the neutron star below which the magnetic energy 
density (i.e. magnetic pressure) dominates the ram pressure in 
the accretion disk ( [Lamb et al.|1973[ l. 



2 



(27) 



with B = Bns ■ (Rns/R)^ the magnetic induction of the neutron 
star (NS) in vacuum, Bns the stellar surface field at the magnetic 
equator, p the average density at radius R and vk = VGMns/^ 
the Keplerian orbital velocity. 

We use the spherical continuity equation (see justification in 
Sect. |3A3] l 



Mns - AnR pvf{. 



(28) 



with Vff = V2vk the free-fall velocity, to eliminate p and solve 
for the Alfv en radius ( [Davidson & Ostriker] [T973) [Eisner &[ 
[Lambl 1977) 1 



2GMnsM2J 



1/7 



where jj - Bns^^s '■^^ magnetic dipole moment of the neutron 
star. 

Since the magnetic field lines are dragged along with the 
neutron star's rotation, their orbital angular frequency equals the 
spin angular frequency Ws of the neutron star. At a distance R 
from the neutron star, the corotating field lines have an orbital 
velocity cOf^R. This equals the Keplerian velocity at the corota- 
tion radius 



Rc, 



/GMns\ 
I ^1 I 



1/3 



(30) 



The propeller effect acts only if the magnetic field energy domi- 
nates the kinetic energy of the disk beyond the corotation radius, 
and when the orbital velocity of the field lines exceeds the es- 
cape velocity at the Alfven radius ( Rappaport et al.[r2004[ ), i.e. 



LOsRfi > V2vk.;,, where vk.^^ is the Kepler velocity at 7?^, so 



(jjs > V2(yK,/j which implies 7?^ > l^^^R^o. Inserting Eqs. (29i 
and (|30| yields the propeller criterion 



M, 



NS 



(2GMns)5/3 ■ 



(31) 



3.4.3. Disk accretion 



By using the thin disk description by Dunkel et al. (2006) it fol- 
lows that the magnetic field is far too weak to dominate the gas 
flow in the disk even for the lowest expected accretion rates. 
Even for the slightly unrealistic case of a neutron star with a 10^ 
G equatorial magnetic field and a 1 ms spin period, the mag- 
netic field only dominates the kinetic energy of the Kepler flow 
at the corotation radius if the mass transfer rate is lower than 
10"'^ Mgyr"', ~ 10'" times lower than for the spherically in- 
falling case. High density blobs of gas can also dominate the 
magnetic field and be accreted (Aly & Kuijpers 1990) . Such 
blobs could form if matter piles up near the Alfven radius, held 
back by the magnetic field, and enter once enough matter has 
accumulated. 



The conclusion is that Eq. ( 29 1 is an upper limit for the 
Alfven radius since the accretion is not spherical as assumed in 
Eq. ( [28] l but occurs in a disk. However, since the density and the 
geometry of the interacting disk-magnetosphere region is poorly 
known, we choose to allow for disk accretion in the equations by 
also considering an effective magnetic field strength at the stellar 
surface which is smaller than the real field strength Bns. 

3.5. Disk Instability Model 

[in't Zand et al. ( [2007| l make a rough estimate for the lowest sta- 
ble mass transfer rate for helium composition based on Lasota 
(2001) and [Menou et"aL] ( |2002| ) 

vO.3 



3 ■ 10 



-10 



\Mo/ I hr / 



Mgyr 



(32) 



where f orb is the orbital period. Once the mass transfer rate de- 
creases below this value, a thermal-viscous instability appears 
(Sect. [2.6.4| l and the disk experiences outbursts. The timescale 
on which an accretion disk rebuilds after a collapse is equal to 
Mdisk/Mii, where the disk mass is given by 



(29) M, 



disk 



r 



2nRI.dR, 



(33) 



with S the surface density for zone C in Dunkel et al. ( 2006| l. Rin 
and Rout are the inner and outer disk radii and will be estimated 
later. 

The Disk Instability Model is still not fully understood, as 
illustrated by the case of SS Cygni, a dwarf nova that accord- 
ing to the Disk Instability Model is too bright and hot to have 
an unstable disk ( [Schreiber & Lasota|2007| l. This increases the 
uncertainty in the mass transfer rate given by Eq. (32i. On the 



other hand, the intermediate state in AM CVn systems seems to 
be predicted quite well ( [Tsugawa &^Osaki| l997| i. 

4. Results 

4. 1 . Feedback of angular momentum 



Since rj is the most uncertain mass flow parameter in Sect. 3.2.1 



here we first look into the likelihood of limited feedback of an- 
gular momentum at some low mass ratio. If the torque between 
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Figures. Important radii within the accretion disk; the dis- 
tance between the accretor and the LI point c/li (solid), the 
average equatorial Roche lobe radius (i.e. in the orbital plane) 
(dashed), the circularization radius Rdi-c (dash-dotted) and /?circ 
corrected for donor influence during free-fall, (dotted, Eq. 
19 1. Furthermore, the pluses (+) indicate the outer tidally stable 
radius from |Papaloizou & Pringle| ( |1977) l. The stars (★) indi- 
cate the largest orbit that does not intersect others, according to 



Paczynski ( 19771 (which lies near the 3:1 resonance). The two 
open circles are outer disk estimates by |Priedhorsky & Verbunt| 



, and the three filled triangles are outer disk radii from our of Priedhorsky & Verbunt 



SPH simulations. 



outer disk and donor is too weak, a dynamical instability may 

resuh (Sect. |2.6.1| ). 

Also mentioned in Sect. 2.6.1 was that l Yungelson et al.| 
( 2006| l and Lasota (2008i have suggested that a resonance be- 
tween the orbital period and the (mean) period of particles in the 
disk may prevent the accretion disk from redistributing angular 
momentum. Also, when q < 0.02, the circularization radius ex- 



Figure 6. Numerical smoothed particle hydrodynamics simula- 
tion of an accretion disk in dynamical equilibrium for mass ratio 
lO^^* after 400 orbits. The solid curve represents the equatorial 
Roche lobes. The accretor has been marked by a small open cir- 
cle. The dashed circles indicate the 3:1 (inner) and 2:1 (outer) 
resonances with the orbital period. 



numerics). The simulations show that the disk fits well inside the 
Roche lobe (triangles in Fi g. [5]l. This confirms the conclusion 
' " ' ' ' ' 1988| l, who used an accretion disk 



ceeds the estimates of the outer radius by Paczynski ( 1977 1 and 
[Papaloizou & Pringle ( 1977) (~ 0.48fl, see the stars and pluses in 
Fig.|5]) which could hinder accretion, and at the same time reduce 
the feedback of angular momentum, making the binary less sta- 
ble. Their value of ^ = 0.02 can be refined by using the adjusted 
circularization radius equation by Verbunt & Rappaport ( 1988| l 
(dotted curve in Fig. E\, who correct the circularization radius 
Rchc/a = (1 +q)(di^i I of ( [Frank et al.|2002| (dash-dotted) for the 
effect of the donor on the stream of matter During this process, 
angular momentum is extracted from the stream, so R\i < R^ 
Using (Eq. T9ll instead of lowers the critical mass ratio to 



q - 0.01 (Fig. 5f7Figure|5]also shows that the location of the first 
Lagrangian point (solid curve) is very close to the donor at low 
mass ratio, and that the equatorial Roche lobe radius (dashed) 
also approaches the LI point, because the equatorial accretor 
Roche lobe becomes more circular rather than teardrop-shaped 
at low mass ratioQ 

We ran numerical smoothed particle hydrodynamics (SPH) 
simulations such as that in Fig.|6]which show that disks function 
properly even at a very low mass ratio. The figure shows a disk 
in dynamical equilibrium after 400 orbits (see Simpson & Wood 
( [T998] l and |Wood, Thomas, & Simpson, (2009) for details of the 



^ Note that the equatorial Roche lobe radius is significantly larger 
than the average radius over the whole volume (Eq. |12[ l, similar to an 
equatorial bulge. 



model to find that the disk expansion is not stopped at 0.48fl 
and that the disk stays within the accretor Roche lobe (circles in 
Fig. |5]) and tidal torques are strong enough to transfer back all 
angular momentum. Also, nothing special appears to happen in 
the disk near the most important resonant radii, the circles in Fig. 
|6] Hence, 77 remains zero at very low mass ratio. In the remainder 
of the paper we will consider the full-feedback case, except for 



Sect. 4.2.3 where we will show evolutionary tracks for the case 
in which feedback somehow stops at a low mass ratio. 

Observations of binaries cannot rule out the existence of a 
dynamical instability since no systems with an upper mass ratio 
below 0.01 have been discovered, although a couple have mini- 



mum mass ratios below this limit ( 


Galloway et al.|2002 Krimm 


let al.|2007t|Altamirano et al.|2010 





4.2. The evolution of ultracompact X-ray binaries 

The key elements in understanding the evolution of UCXBs are 
represented by Eq. (j6]l. Mass transfer is driven by gravitational 
wave radiation, and also depends on the relative responses of 
donor and Roche lobe size to mass transfer. In Sect. 14.11 we 
found that advected angular momentum will be returned from 
disk to orbit for any mass ratio, which constrains and allows 
us to solve all relevant binary parameters during the entire evo- 
lution when assuming stable accretion, rather than in outbursts. 
A system can survive the onset of mass transfer only if the ac- 
cretor can eject all transferred matter exceeding the Eddington 
limit via isotropic re-emission. For any realistic stellar-mass ac- 
cretor, a merger can be avoided if the white dwarf donor has a 
mass below ~ 0.3 Mq. Mass transfer is dynamically stable if, 
upon mass loss, the donor expands less than its Roche lobe (Eq. 
|8]l, that is, when the solid curve in Fig. |7] lies below the dashed 
curve. UCXBs with black hole accretors can theoretically ex- 
ist, but may be rarer than neutron star UCXBs due to the steep 
slope of the initial mass function. However, the binary survival 
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Figure?. Mass-radius relation exponent (Sect. 3.2 1 for var- 
ious modes of mass transfer from a white dwarf to a neutron 
star of initially 1.4 with a 12 km radius. The lower dotted 
curve shows the case of accretion efficiency e = and full feed- 
back (77 = 0) (almost all transferred matter lost via isotropic re- 
emission, hardly any accretion). The middle dotted curve is for 
e = 1 and 77 = (conservative mass transfer). The upper dot- 
ted curve is for no feedback (77 = 1). The zero-temperature white 



dwarf exponent (Sect. 3.3 1 is represented by the dashed curve 

0.59 (Ma 



23| l 
30th 



The soHd curves starting at ^ = 0.59 {Ma - 0.83 Mq, Sect, 
indicate the actual value of (for e shown in Fig. |4|i for 
the full feedback and no-feedback case. Until q - 0.01 there is 
full feedback in either case. Below q - 0.01 the full feedback 
case is shown by the lower solid branch, where the open disk 
denotes an age of 10 Gyr (since the onset of mass transfer). The 
no-feedback case has a jump in at = 0.01, after which it fol- 
lows the upper dotted curve until the intersection with the dashed 
curve (filled disk) where donor disruption occurs at a system age 
of 300 Myr. 



rate during the supernova event can be different for black hole 
formation than it is for neutron star formation due to a different 
kick velocity distribution. 

4.2.1 . Evolution in tine case of full feedback 

In Fig. [8] the mass transfer rate starts very high because of 
the massive donor and especially because of the short orbit. 
The first part of the tracks lies above the isotropic re-emission 
limit (dashed, Sect. |2.4| i, which means that a system cannot sur- 
vive here, this part of the tracks is shown just for illustration. 
Isotropic re-emission of super-Eddington mass transfer is pos- 
sible for donor mass below 0.37 and orbital period above 
1.8 min. For lower donor mass (0.08 Mq, with a corresponding 
orbital period of 8 min), the mass transfer rate decreases below 
the Eddington limit (dotted), and from here on conservative mass 
transfer is possible (Fig. |4|i. At later times, the mass transfer rate 
reaches the thermal-viscous disk instability rate (dash-dotted. 
Sect. |3.5| l. The donor mass is about 0.02 Mq at this point, and 
the orbital period 28 min in the case of a helium disk. Systems 
with a longer period than this do no longer have stable accretion 
disks. After 10 Gyr of mass transfer, a system with a 1 .4 M© neu- 
tron star accretor would have a donor mass of 4.2 ■ 10"^ Mq and 
an orbital period of 85 min. In the case of a 10 Mq black hole 
accretor, the donor mass would be 2.6 ■ 10"^ Mq and the orbital 
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Figures. Mass transfer rates for a helium white dwarf donor 
with an initially 1.4 Mq, 12 km radius neutron star accretor 
(lower solid curve) and a 10 Mq black hole accretor (upper solid 
curve), starting at the dynamical instability limit (i.e. ignoring 
possible merging), for full feedback. The dashed curves branch- 
ing away from these represent the evolution assuming no feed- 
back below q - 0.01 (Sect. |4T| . The numbers correspond to the 
two circles on their left-hand side and represent the logarithm of 
the system age (yr), which is the time since the onset of mass 
transfer. Indicated are the isotropic re-emission regime (dashed, 
Eq. Ill, the Eddington limit (dotted) and the accretion disk in- 
stabiHty limit by i n't Zand et al.| ( p007 ) (dash-dotted. Sect. [33 1 
Note that the dashed neutron star curve at its local maximum 
'touching' the upper dash-dotted curve is a coincidence. 



period 111 min. (The 10 Gyr given here, as well as ages given 
elsewhere in the paper, is the time measured after the onset of 
mass transfer to the neutron star or black hole.) These systems 
technically do not fall under the present (observational) defini- 
tion of UCXBs (an orbital period < 60 min). Gravitational wave 
radiation is stronger in systems with a more massive accretor, 
therefore the black hole track in Fig. [8] lies above the neutron 
star track (to be explained in Sect. |4.2.2| i. 

Figure |9] is a linear graph showing more clearly how the or- 
bital period increases with time, initially the period increases 
very rapidly, and slowly later on. Systems with a massive ac- 
cretor reach a given orbital period earlier on, because of their 
shorter evolutionary timescale. 

The change in accretor mass during evolution For young sys- 
tems with massive donors, the orbit is so compact that the donor 
is forced to lose mass at a high rate in order to avoid a merger. 
For these systems e 0, because most mass is expelled from 
the system by isotropic re-emission. Once the transfer rate drops 
below the Eddington limit, e- \ from Eq. ( 17 1 as all transferred 
mass can be accreted (Fig.|4]i. This happens at a donor mass of 
0.08 Mq. Because e goes from negligible to 1 quickly (in terms 
of donor mass), the total mass increase of a neutron star AMa 



will not exceed - J eMd df = 0.08 Mq by Eq. ( 16 1. The angular 



momentum associated with AMa is relevant when considering 



spin up of the accretor, in Sect. 4.3 



Observed potential long-period UCXBs The binary millisec- 
ond pulsar HETE J 1900. 1-2455 has an orbital period of 83.3 
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Figure 9. Orbital period against time for a helium white dwarf 
donor with an initially 1.4 Mq, 12 km radius neutron star (lower 
solid) and 10 Mq black hole (upper solid) accretor. 



min, but a donor much more massive than systems of the same 
period following the evolution described here. The donor is 
probably a brown dwarf of 0.016 - 0.07 Mpl lKaaret et al.|2006| l. 

Because of their shorter evolutionary timescale, white dwarf 
- black hole systems can potentially reach an orbital period of 
1 10 min within the age of the Universe, if the system is not dis- 
rupted at some low mass ratio and mass transfer starts relatively 
quickly after the Big Bang. There are 2 low-mass X-ray bina- 
ries known with 90 - 115 min periods, viz. IE 1603.6+2600 
at 111 min (Morris et al.|[T990l l and XTE 1748-361 at 97 min 



( Bhattacharyya et al. 2006[ l. Another 2 low-mass X-ray bina- 



ries have periods just outside this range, SAX J1808.4— 3658 at 
120 min (|Chakrabarty & Morgan |[i^8| and GS 1826-238 at 
126 min ([Homer et al.|1998| l, although the last period is uncer- 
tain ( |Liu et al.||2007| l. All these systems exhibit type I X-ray 
bursts so they must harbor a neutron star. The donors are not 
expected to be cold white dwarfs since these long periods can 
not have been reached by evolution on the gravitational wave 
timescale. 



Why is conservative mass transfer less stable than non- 
conservative mass transfer? for conservative mass transfer 
(e = 1) is higher than for fully non-conservative mass transfer 
(e = 0) (the lower two dotted curves in Fig. |7]i implying dy- 
namically less stable mass transfer for e = 1 by Eq. (j6]l, when 
we assume isotropic re-emission and full feedback of angular 
momentum from disk to orbit (77 = 0). This may seem counter- 
intuitive since less mass ejection makes a system more stable 
because less angular momentum is lost. The explanation is that 
the latter effect (represented by the y-term in (a, Eq. (|24]l) is 
outweighed by the effect of the (1 - /3)q (= eq here) term for 
all e and q due to the unfavorable effect on the accretor mass: 
in the case of non-conservative mass transfer the accretor gains 
less mass, implying a must increase more to compensate for this 
(7^. = GflM^M^/Mtot) than in the conservative case. This effect 



is enhanced by the increase of with e (Eq. 14 1. For e = 1, 
the mass ratio q = M^/M.^ decreases faster thanTor e = 0. The 
more extreme q causes the donor Roche lobe to shrink more (or 
expand less) upon mass loss relative to a. 



4.2.2. Why do massive accretors speed up evolution? 

Figures [T] and [8] show that UCXBs evolve faster if they have a 
more massive accretor. This is because evolution is driven by 
gravitational wave radiation. Here we give the explanation and 
look into some consequences. 



The effect of accretor mass on mass transfer rate: 
For q < 1, the orbital period of a Roche-lobe filling 
system is only weakly dependent on the accretor mass 
(which can be see n from the Roche-lobe approximation by 

for 



fl/3 



Paczyriski ( 1971 1), so the semi-major axis a oc M^^^ 
hxed donor mass, according to Kepler's third law. F rom Eq. 
^ and (7/7)r.wR = -(32/5)(G^ /c^)MdM^M,ot/a'^ ( |Landau 
(^Lifshitz||l975 l follows that Md oc M^ImH^, so the mass 



transfer rate mcreases with accretor mass for a fixed donor, 
even though the orbital separation increases as well. This ex- 
plains why the black hole track in Fig.|8]lies above the neu- 
tron star track. The correction from - is small unless 
(a - approaches zero, as evidenced by the intersection of 
the tracks in Fig. |8] where the black hole system is still fur- 
ther away from its dynamical stability limit (1.04 Mq). In 
reality, this intersection is not important since systems with 
an orbital period below ~ 2 min cannot survive. 
Systems of the same age after the onset of mass transfer: 
When we compare systems of the same age, mass transfer is 
lower in the case of a more massive accretor (again, except 
when they are extremely young). The influence of the lower 
donor mass outweighs the influence of the more massive ac- 
cretor in the gravitational wave equation, causing a lower 
mass transfer rate for more massive accretors of a fixed age. 
An old 1 .4 Mq neutron star accretor system has a mass trans- 
fer rate that is ~ 50% higher than a 10 Mq black hole accretor 
system of the same age, as shown by the fits in Appendix 
[a] (If the systems did not start with identical donors, the 
more massive donor will lose mass fast enough to (almost) 
catch up with the lower mass donor, on a timescale much 
shorter than the typical age of the systems.) A system with a 
lower accretor mass consistently experiences a higher mass 
transfer rate, but still its donor remains more massive. It can 
never catch up, because its mass transfer rate would drop be- 
low that of the black hole system in the case of nearly equal 
donor masses. Lastly, the gravitational wave timescale for a 
system with a neutron star accretor is ~ 7% longer than that 
of a 10 Mq accretor system of the same age, and the orbital 
period is ~ 25% shorter. 

To summarize, in the case of an identical donor, |M(j/M(j| 
and |M(j| are higher for systems with a higher accretor 
mass. For systems of the same age, however, Md and |Md| 
are lower for a higher accretor mass, whereas |Md/Md| is 
marginally higher (i.e. mass is transferred on a slightly 
shorter timescale). 

Consequences for AM CVn systems: 

The above also implies that UCXBs evolve on a shorter 
timescale than AM CVn systems (which have white dwarfs 
both as donor and accretor) with an identical donor The 
higher value of e in AM CVns makes the mass transfer dy- 
namically less stable and the mass transfer rate higher, but 
the accretor mass has a larger impact, unless the donor mass 
is close to the dynamical stability limit. For a low donor 
mass, ^d - is almost independent on the mass ratio and e, 
so the mass transfer rate only depends on the accretor mass. 
Also, AM CVns are expected to experience a higher mass 
transfer rate than UCXBs of the same age after onset of mass 
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Figure 10. Stability regimes in terms of zero-temperature helium 
white dwarf donor mass versus accretor mass. The upper solid 
curve is the upper limit for dynamical stability in the presence 
of an accretion disk that feeds back angular momentum, where 
almost all transferred matter is lost via isotropic re-emission. 
The upper dashed curve is the upper stability limit for isotropic 
re-emission and the lower dashed curve indicates the proposed 
q = 0.01 lower mass ratio limit for feedback (Sect. 4.1 1. The 
lower solid curve encloses the parameter space where systems 
are dynamically stable in the absence of angular momentum 
feedback; all matter enters and stays in the disk. The dash-dotted 
lines indicate when the propeller effect comes into play accord- 
ing to Eq. ( [3T] i for several neutron star spin periods and an equa- 
torial magnetic field of 10^'^ G. The dotted isochrones indicate 
the donor masses for system ages (since the onset of mass trans- 
fer) of 0.1, 1, 10 Gyr based on full feedback. The dynamical sta- 
bility and isotropic re-emission rates are independent of accretor 
radius and type, hence these curves are continuous. Moreover, 
the dynamical stability rates are independent of accretor type 
and therefore also hold for white dwarf accretors. 



transfer, except again for very young systems. Old AM C Vns 
have a more massive donor than equally old UCXBs. 



Dynamically stable parameter ranges Figure [T0| shows which 
combinations of donor and accretor masses are dynamically sta- 
ble, both in the case of feedback and the case of no feedback. 
The accretor mass has a large influence on when an instability 
occurs, if at all. The bottom left region enclosed by the solid line 
is dynamically stable for mass transfer, even in the absence of 
feedback. Initial mass transfer is dynamically unstable for every 
Md if Ma > 1 1 Mq, because of high specific angular momentum 
of the stream in low mass ratio systems. The upper border of this 
region shows that the highest donor mass that can possibly avoid 
dynamical instability in the case of no feedback is 0.42 Mq, for 
an accretor mass of 5.3 MqJ^As explained in Sect. 



2.5 



m prac- 
tice the upper border will be replaced by a limit of ~ 0.5 Mq. 
Note that the donor mass at which the low mass ratio dynami- 
cal instability could occur is very sensitive to accretor mass near 
Ma = 1.4 Mq. At exactly which donor mass this instability oc- 
curs depends on the mass-radius relation of the donor, which in 
turn depends on its temperature and composition. 

The upper limit for dynamical stability in the case of full 



10 



feedback is represented by the upper solid curve in Fig 
Because (Eq. \T3\ has a minimum of -5/3 for q ^ regard 
less of e and y, any system with a zero-temperature white dwarf 
donor mass exceeding 1 .07 Mj^^is dynamically unstable, even for 
arbitrarily high accretor mass 



Age of disruption If feedback would suddenly stop at ^ = 0.01, 
a dynamical instability would happen immediately if Ma > 
4.7 Mq, and later for less massive accretors (shown by the inter- 
section between the lower solid and dashed curves in Fig 



lOi 



For an UCXB with a 1 .4 Mq accretor, the instability would hap- 
pen at an age of 300 Myr, for Md = 3.2 ■ 10"^ Mq (the upper-left 
intersection at ^ = 2.3 ■ lO^^* in Fig.|7]i. The reason such a low 
donor mass could be reached well within 10 Gyr is the increased 
mass transfer rate caused by the extra angular momentum sink, 
about one order of magnitude for a 1.4 Mq accretor (Fig. [8]). A 
10 Mq black hole system would reach q - 0.01 after already 80 
kyr. 

Disruption of low-mass donors would happen later than 
10 Gyr if Ma < 0.6 Mq, this could apply to AM CVn systems. 
UCXBs with any realistic accretor mass would suffer disruption 
well within 10 Gyr. 



4.2.3. Evolution in the case of no feedback at low mass ratio ^'^^ ^^^f^'" ^^^^^ P'°P^"^' P^^'^ 



In Sect. 2.6.1 we mentioned several proposed mechanisms for 
reducing feedback of angular momentum from the disk to the 
orbit, but we think that none of these will actually work (Sect. 
|4.1| l. In the case that after all a mechanism exists that leads to 
a reduction in or absence of feedback in systems below a cer- 
tain low mass ratio, here we briefly discuss the evolution of such 
systems, including the dynamical instability that may follow. We 
consider the case in which all feedback stops once q < 0.01, i.e. 
T] = 1 and j0 = in Eqs. (24 1 and ( 15 i for these mass ratios. 



Under this condition, the evolution of ^ for a neutron star ac- 
cretor UCXB is shown by the upper solid curve in Fig. |7] The 
sudden jump in ^ leads to an increased mass transfer rate per 
Eq. (|6]), which speeds up evolution, therefore the low mass ra- 
tio where ^d and meet can be reached well within the age of 
the Universe. Two evolutionary tracks are shown as the dashed 
curves branching away at <7 = 0.01 from the solid curves in Fig. 

m 



Even though most likely no dynamical instability occurs at low 
mass ratios, the regular accretion rates as described in Sect. 4.2. 1 
are idealized. Here we investigate the UCXB evolution taking 
a magnetized neutron star accretor into account, using angular 
momentum and energy considerations. This gives insight into 
the feasibility of the propeller effect (Sect. 3.4 1 removing matter 
from the system and hence suppressing accretion onto and X-ray 
emission from UCXBs. 



This upper limit is a result of the function in Eq. 1 13 i with /3 = 
e = having a minimum of = -0.48 at g = 0.080, and (a = -0.48 is 
solved by = 0.42 Mq in Eq. ( (26} . 

' Around an intermediate mass or supermassive black hole, an accre- 
tion disk can only exist if M.^ < 10^ Mq in the case of a 0.6 Mq white 
dwarf donor and Ma < 10* Mq in the case of a 0.1 M© white dwarf 
donor, since otherwise the innermost stable circular orbit exceeds the 
semi-major axis. 
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Figure 11. Radii in UCXB accretion disk versus orbital period 
for a 1 .4 Mq neutron star accretor with a residual equatorial mag- 
netic field of 10** G. The shaded layer at the bottom represents 
the neutron star surface, at a radius of 12 km. The curve an- 
notated LI shows the distance of the first Lagrangian point to 
the neutron star, which is an upper limit on the outer disk radius. 
The circularization radius by |Verbunt & Rappaport| ( [T988j ) is dis- 



cussed in Sect. 4.1 The two upper dashed curves indicate the 
radii where the temperature of a helium disk (Dunkel et al. 2006| l 
equals the He I ionization temperature, for two values of the vis- 
cosity parameter a, as discussed in Sect. 



2.6.4 The solid coro- 



tation radius curve applies to the full propeller case, the dashed 
one to the hypothetical full accretion case. The speed-of-light 
cyUnder radius and Alfven radius are also shown. 



4.3.1 . Steady mass transfer 



First we assume the absence of disk instability (Sect. 3.5 i. This 
means that matter approaches the magnetosphere at a constant 
rate (that is, only varying on the evolutionary timescale). From a 
magnetic perspective, the evolution can be divided in a spin up 
and a spin down stage. 

Spin up Initially the mass transfer rate is so high that the ki- 
netic energy in the disk completely dominates the magnetic field 
energy; the Alfven radius lies either at the surface or just outside 
the neutron star. Moreover, the relatively low initial accretor spin 
frequency (and correspondingly large corotation radius) ensures 
that any interaction between disk and magnetic field that is still 
present, causes the neutron star to spin up rather than matter to 
be expelled, because the fast rotating disk applies a torque on the 
field lines, which leads to transfer of angular momentum from 
disk to neutron star. We assume that as long as the Alfven radius 
is smaller than 2'^-^ ^ 1.26 times the corotation radius (Sect. 
|3.4.2| l, all arriving matter is accreted, and the amount of angu- 
lar momentum added to the accretor corresponds to the specific 
angular momentum of disk matter at the Alfven radius. Specific 
angular momentum of matter in a Keplerian orbit is proportional 
to the square root of the distance to the accretor, and we will 
find that the typical ratio of the outer disk radius to the Alfven 
radius exceeds lO'' during the spin-up stage. Therefore, the frac- 
tion of angular momentum that is not transferred back via the 
outer disk is at most a few percent, and our initial assumption of 
dynamically stable mass transfer holds. 



As the neutron star spins up, the corotation radius approaches 
the neutron star At the same time the Alfven radius moves out- 
wards as a result of decreasing mass transfer. At an orbital period 
of 10 min (and a system age of 1 Myr) these two meet. The sur- 
face magnetic field strength at this time is determined by how 



much mass the accretor has gained up to this point (Sect. 4.2.1 
We adopt the empirical relation for the accretion-induced field 
decay by .Shibazaki et al. ^989 ), with a residual field strength 
of Bns = 'lO°-^ G |Zhang & Kojima 2006). We also assume that 
the neutron star does not spin when the mass transfer starts, i.e. 
Jns = 0, since the earlier common-envelope stage is too short- 
lived to allow for significant accretion. Also, the strong initial 
magnetic field suggests a very low initial spin frequency relative 
to the eventual spin frequency. 

The amount of angular momentum added to the neutron star 
by accretion is 



d^NS = dMNs ^JGM^sRu 



(34) 



where Rm is the inner disk radius, which determines the specific 
angular momentum of matter in a Keplerian orbit at this radius. 
We assume that matter stops losing angular momentum outwards 
once it moves inside this radius. Ri„ - Rf^ ( Li & Wang| l999|l ex- 
cept when the magnetic field does not dominate the Kepler flow 
anywhere, in which case the stellar surface is taken, or when Rf^ 
lies outside the speed-of-light cylinder radius Rc - clu^, where 
corotating field lines would travel at the speed of light. There, 
field lines are open and particles can escape by following them 
outwards. The value of the inner radius is summarized as 



/?i„ = min(/?c,max(/?^,/?Ns))- 



(35) 



From 7ns the spin angular frequency Ws is found using /ns = 
/nsWs, where /ns ~ 0.357Mns^n s the m oment of inertia of a 
1.4 Mq, 12 km neutron star ( Lattimer & Schutz 2005 ). We con- 
sider angular momentum accretion rather than energy accretion, 
because the latter is complicated by radiation losses and orbit 
coupling. 



Spin down Once the Alfven radius moves beyond the corota- 
tion radius, the torque is reversed, and matter in Keplerian orbits 
between the corotation radius and the Alfven radius may be ac- 
celerated by the super- Keplerian velocity of the field lines, re- 
ducing the accretion rate. Matter near R^ (or R^ once smaller) 
that reaches the escape velocity then leaves the system on tra- 
jectories outside the orbital plane, the angle between the out- 
flow and the plane being smaUer for a stronger propeller effect 
( [Lovelace et al.|1999t[Ustyugova et al.|2006| l. 

Due to its accretion history, the neutron star has gained a 
large amount of rotational kinetic energy at this stage. The evolu- 
tion of spin frequency u)^ once R^ > l^^^R^o depends on whether 
there is enough energy, and whether it can be employed to un- 
bind all or some fraction of arriving matter. At this point we 
assume that there is sufficient energy to unbind all transferred 
matter when Rf, > l^^^R^o, and will later find that this is indeed 
the case. Then, the energy of the neutron star changes as 



d^NS 



2R\„ 



dt. 



(36) 



where d^NS and are negative, and the factor 2 appears be- 
cause the total energy in a Keplerian orbit is half the poten- 
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tial energy. The rotational kinetic energy of the neutron star is 
-Ens = j/nsw^, which yields 



Spin evolution Figure 11 illustrates the evolution of the accre- 
tion disk, magnetosphere and corotation radius for an UCXB. 
The dotted part of the Alfven radius curve (Eq. [29) 1 lies beyond 
the speed-of-light cylinder radius, and therefore the field lines 
cannot corotate in this region. The corotation radius is given by 
Eq. (30 1 and is proportional to p1^^ , where the accretor spin pe- 
riod f s = 27T/a>s- The solid corotation curve shows the evolution 
when no matter is accreted once > 2''^7?co, but instead all 
transferred matter is unbound from radius Ri„, extracting energy 
from the rotation of the neutron star. For comparison, the dashed 
corotation curve shows the spin evolution in the hypothetical 
case in which, once > l^^^R^o, aU matter would continue to be 
accreted with specific angular momentum corresponding to Ri„ 
instead of stopped by the magnetosphere. Before - l^^^R^o, 
both assume full accretion so there these two are equal. 

Because the corotation radius stays within the Alfven radius, 
apparently the early accretion at a high rate has sufficiently spun 
up the neutron star to provide enough energy to unbind any frac- 
tion of transferred matter during the whole remainder of the evo- 
lution, when the mass transfer rate is low. This justifies our as- 



sumption above Eq. ( 36 1 



By using a lower effective equatorial surface magnetic field 
strength to accommodate for disk accretion, as discussed in Sect. 
3.4.3 less angular momentum is accreted per unit mass due to 
the smaller Alfven radius, but also accretion can continue for a 
longer time because the propeller effect appears at a later stage. 
The result is that the neutron star is spun up to a lower spin pe- 
riod, while subsequently less energy is needed to unbind arriving 
matter because the propeller effect starts at a lower mass trans- 
fer rate. The eventual spin period is lower in the case of a lower 
effective magnetic field. 



Figure 12 shows the spin period directly, rather than the re- 
lated corotation radius. The hypothetical case in which the pro- 
peller effect does not work, so all matter is accreted at the stellar 
surface carrying the local Keplerian specific angular momentum, 
is shown by the dash-dotted curve (Sect. 3.4.3 1. For comparison 
we include the observed UCXB neutron star spin periods, which 
are too small in number and have too scattered spin periods to 
have any constraining value to our model. 

Now we know that full propeller is energetically possible for 
UCXBs, does it actually happen? The observed neutron star spin 
periods, typically 5 ms, suggest that at least some spin down 
is plausible. On the other hand, the faint transient X-ray binary 
SAX J 1808.4-3658 proves that at least some accretion is still 
possible onto a very fast spinning (2.5 ms) neutron star, since 
it shows typ e I X-ray bursts (|in 't Zand et al.|[T998l ). [Spruit & 



Taam 



( |1993| l; [Rappaport et al. 



( 2004 1 showed that some accre- 



tion likely persists assuming the inner disk radius comes quite 
close to the corotation radius. They also mentioned that this disk 
structure probably breaks down at very low mass transfer rate, 
and that full propelling probably happens then. Simulations by 
D'Angelo & Spruit (2011 1 showed that accelerated matter does 



not necessarily leave the system, because the inner disk radius 
may remain trapped at some short distance outside the corotation 
radius. This scenario, however, appears more likely for stronger 
magnetic fields (~ 10'^ G) than the ones considered here, and 
for stars with a much shorter spin-down timescale than millisec- 



The solution for 0)^ is implicit since a)^^ depends on £ns, which itself 
depends on the speed-of-light cylinder radius via Eqs. (|36b and (|35|. 
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Figure 12. Spin period versus orbital period for a 1 .4 Mq neutron 
star accretor with a residual equatorial magnetic field of 10^^ G, 
assuming all transferred matter is pro pelled out of the system 



once Rf, > l^'^R^. 



(solid. Sect. 



3.4.2 1, or hypothetic ally keeps 



being accreted (dashed), or hypothetically is accreted at the neu 
tron star surface all the time (dash-dotted). Observed millisecond 
pulsar-UCXB parameters are represented by filled disks; XTE 



J 1807-294 at an orbital period of 40.07 min ( Markwardt et al 
2003a) and a spin period of 5.25 ms (Markwa rdt et al.|2003b| l, 



XTE J175 1-305 at 42.42 min, 2.30 ms (Markwardt et al. 



200211, 



XTE J0929-314 at 43.58 min, 5.40 ms ( [Galloway et al.|20021i, 
SWIFT J1756.9-2508 at 54.70 min, 5.49 ms (Krim m et al.[2007ll 
and NGC 6440 X-2 at 57.3 min, 4.86 ms (Altamirano et al.[ 
2010). The open disk represents 4U 0614-1-091, which harbors 
a 2.41 ms pulsar ( jStrohmayer et al. 2008[l and has a tentative 
orbital period of 51.3 min ( [Shahbaz et aL[[2008| [Hakala et al] 
|20TT1 ). 



ond pulsars. Therefore, a non-accreting state seems more likely 
for UCXBs containing a recycled neutron star, and the propeller 
effect may work depending on whether the inner disk radius ex- 
ceeds the corotation radius by a sufficiently large ratio. 

Whether accretion at very low mass transfer rate is possible 
or not in a system with a stable accretion disk perhaps is not the 
right question to ask, since accretion most likely occurs irregu- 
larly. Figure [TT suggests that a helium disk (and also a carbon- 
oxygen disk, which is not shown) remains divided in a neutral 
outer part and an ionized inner part, where the inner part at all 
times is too large for the magnetosphere to disrupt. This implies 
that the disk instability remains relevant and will come into play 
once the mass transfer rate becomes too low (Sect. [2A4] l. 

4.3.2. Non-steady mass transfer 

Varying mass flow rates and densities near the accretor have a 
large impact on the interaction between the disk and the mag- 
netic field. For instance. Spruit & Taam ( 1993 i described cycles 
of accretion, during which matter piles up outside the Alfven ra- 
dius until the local density is such that the inner disk forces the 
Alfven radius inwards, allowing accretion. When the disk den- 
sity decreases, the Alfven radius moves outwards again. In Fig. 
1 1 it can be seen that the Alfven radius at Eddington-limit mass 
transfer is almost equal to the late-time corotation radius, but 
taking into account the decreasing magnetic field, the Alfven ra- 
dius will always lie well inside the corotation radius at the maxi- 



14 



van Haaften et al.: The evolution of ultracompact X-ray binaries 



3.0 
2.5 

g 2.0 
M 1.5 



S 1-0 



o 0.0 
-0.5 



-1.0 





NS, a = Ci.Vl^^^ 




^^^^3^BH, a = 0.02 






- 


^ , ->&^^NS, ff = 0.1 ; 




' ^^^^ - 




BH, ff = 0.1 


/ / 
/ / 









20 



40 60 80 

Orbital period (min) 



100 



120 



Figure 13. Recurrence time for a 1 .4 Mq neutron star accretor 
with a residual equatorial magnetic field of 10^'^ G (solid) and 
for a 10 Mq black hole accretor (dashed) with a helium disk 
described by Dunkel et al. (2006), for two values of the viscos- 
ity parameter a, of which 0.02 probably best describes the disk 
viscosity in quiescence. Disk instability becomes relevant to the 
right of the vertical line at 28 min. The circles indicate an age of 
10 Gyr The discontinuity in the slopes of the neutron star curves 
at an orbital period of 29 min is because the speed-of -light cylin- 
der defines the inner disk radius beyond this period (Fig. [TT] ). 



mum accretion rate, so a near-Eddington outburst can always be 
accreted. 

In Sect. 13.51 we have discussed under which circumstances 
accretion could become irregular due to thermal-viscous disk in- 
stability. When the mass transfer rate has become sufficiently 
low, the disk periodically collapses and builds, which leads to 
high mass transfer rates near the accretor every cycle. The re- 
currence time is the time it takes until the empty disk is refilled. 
Figure [13] shows the recurrence time for our UCXB model. The 
inner disk radius is given by Eq. ( 35 1 in the case of a neutron 
star accretor and 3 Schwarzschild radii in the case of a black hole 
accretor The outer radius is estimated by R^vx ~ 0.9iiLi based on 
Fig.|5] The recurrence time is increasing because the mass trans- 
fer rate decreases faster than the disk mass. During quiescence 
the disk viscosity is low, of the order of o- = 0.02 ( |Lasota|2001| l. 
Systems of age 10''"'" yr with a neutron star accretor have a re- 
currence time of hundreds of years. In the case of an equally old 
system with black hole accretor, this time is a bit shorter The 
long recurrence time suggests that systems are only rarely ob- 
servable as bright X-ray sources, so that at any given time, only 
a small fraction of a population is visible. 



5. Discussion and conclusions 

We have studied the evolution of ultracompact X-ray binaries, 
and in particular their characteristics at an old age, when they 
reach a very low mass ratio and low mass transfer rate. The aim 
is to investigate how likely the occurrence of a dynamical insta- 
bility is, and how important the propeller effect and disk insta- 
bility are. 

If the feedback of angular momentum from disk to orbit were 
to be reduced, the lifetime of an UCXB would be dramatically 
shortened. Based on accretion disk simulations (Sect. |4.1| i, we 
conclude that an accretion disk in a system of any realistically 



achievable mass ratio is neither significantly affected by reso- 
nant orbits, nor does it expand outside the accretor Roche lobe. 
The latter is in agreement with Priedhorsky & Verbunt ( 1988) . 
Most likely the disk is still capable of full feedback of angular 
momentum to the orbit so that the donor will never be disrupted. 
Consequently, unless a new mechanism is discovered that limits 
or stops feedback, or pulsar radiation evaporates the donor by 
heating and stripping off the outer layers (Kluzniak et al.|19"88 



IRuderman efa l. 1989; Shaham|1992| i, or UCXBs become tem 
porarily detached at some stage, the vast majority of UCXBs 
today will have orbital periods around 70 - 80 min and experi- 
ence low mass transfer rates of 10"'^ - 10"'^ Moyr"' , depending 
on the precise mass-radius relation of the donor (Systems that 
started mass transfer more recently will have shorter orbital pe- 
riods, but since the evolutionary timescale increases rapidly with 
increasing orbital period, most systems are expected to have a 
long orbital period at the present. The time between the zero-age 
main sequence and the onset of mass transfer is in general short, 
as will be shown in a forthcoming paper.) 

However, all ~ 15 known UCXBs have much shorter or- 
bital periods, below 60 min. Since these short-period systems 
evolve on a much shorter timescale than long-period systems, 
there must exist a much larger population of apparently invisi- 
ble, old, long-period UCXBs in the Galaxy. A strong selection 
effect hides the large population of relatively common systems 
so that we see only a couple of much rarer systems during their 
short-lived high mass transfer phase. The bolometric luminosity 
corresponding to a constant mass transfer rate of ~ 10"'^^ Moyr"' 
is ~ 5 ■ lO-'^ erg s"', within the observational limits of the best 
present-day X-ray telescopes even at a distance of 8 kpc (the 
distance to the Galactic Bulge) if these systems emit predomi- 
nantly in X-rays. Then, the only consistent resolutions are that 
these long-period systems are even fainter than expected, or vis- 
ible only a small fraction of the time. These resolutions are sup- 
ported by the propeller effect and the disk instability model, re- 
spectively. 

The propeller effect probably strongly influences accretion in 
UCXBs, especially at very low mass transfer rates. In the event 
of a stable disk, rotational kinetic energy stored in the neutron 
star during the accreting stage is sufficient to expel all matter 
arriving at the magnetosphere from the system once the propeller 
effect starts to work. Therefore, systems with a stable accretion 
disk could have a significantly reduced X-ray luminosity. 

However, the disk instability model, which describes how 
disks periodically build and collapse, still applies at the very low 
mass transfer rates expected in UCXBs. During a collapse, mat- 
ter is capable of penetrating the magnetosphere and reach the 
stellar surface, which means that observationally the system ap- 
pears as a transient X-ray source. Based on estimates of the disk 
structure, the recurrence time of old (10 Gyr) systems is of the 
order of a few 100 yr, depending on at what mass transfer rate 
the disk becomes unstable, among other things. Since black hole 
accretor systems do not suffer from the propeller effect, these 
would be expected to be better visible, but none have been dis- 
covered so far Apparently they are much rarer than their neutron 
star accretor counterparts, or disk instability is the dominant pro- 
cess in making low mass ratio UCXBs nearly invisible. 

A third mechanism that can inhibit the X-ray visibility of a 
black hole UCXB is the inefficiently radiating disk (ADAF - ad- 
vection dominated accretion flow) ( Narayan & Yi|1994 



Lasota 



2008 ). In ADAFs, most of the thermal energy stored in the inner 
disk is accreted into the event horizon. For a neutron star UCXB 
this does not apply as the energy would still be released upon 
accretion. 
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Lastly, UCXBs with a low mass transfer rate may also be- 
come fainter in X-rays if they emit a significant fraction of their 
energy as UV or even optical radiation. The Dunkel et al. ( 2006| l 
central disk temperature for helium composition is given by 

oc M^^'" or M^l^ depending on the choice of opacity, while 
a thick disk has T oc M^l'^, so the disk spectrum shifts to longer 
wavelengths as the mass transfer rate decreases. The hot inner 
part of the disk, origin of most radiation, still emits mainly X- 
rays, with the result that the disk as a whole emits more than half 
of its energy as X-ray radiation, even at the lowest mass transfer 
rates we find. The inner disk, however, is affected by the mag- 
netic field in the case of neutron star accretors (in particular at 
low mass transfer rate), and can become an ADAF in the case 
of black hole accretors. In these cases the outer disk spectrum 
becomes more important and the system may emit most radia- 
tion in the UV, which is difficult to observe due to absorption 
(jCampana & Stella|2000j). 

The consequences of the evolutionary behavior of UCXBs 
for the population in the Galactic Bulge will be the topic of a 
forthcoming paper 
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Appendix A: Fitted tracks 

Sometimes it is useful to have simple analytic estimates of the 
evolution of the quantities involved. For this purpose we provide 
powerlaw approximations of the tracks with full feedback of an- 
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Table A.l. Values for the parameters in Eqs. ( A.l 





Fitted tracks 


Analytic approximations 


Parameter 


1.4 Mo 


10 Mo 


1.4 Mg 10 Mo 


/i 


1.18 


0.882 


0.736 0.515 


h 


0.326 


0.482 


0.196 0.281 


h 


5.92 


4.62 


2.70 1.89 


/4 


3.01 • 10-3 


1.61 • 10-2 


3.71-10-" 1.38-10-3 


/s 


1.38 


1.69 


0.530 


k 


7.75 • lO"'' 


1.87 • 10-3 


7.15-10-3 2.65 -10-- 


ei 


-1.29 


-1.29 


-14/11 (« -1.27) 


ei 


0.242 


0.237 


3/11 (~ 0.273) 


es 


-0.313 


-0.323 


-3/11 {- -0.273) 


en 


-5.32 


-5.45 


-14/3 {- -4.67) 


es 


-1.29 


-1.36 


-1 


66 


4.11 


4.01 


14/3 (~ 4.67) 



Notes. 1 .4 Mq and 10 Mo are the accretor masses. The fitted tracks are 
discussed in Appendix [A| and the analytic approximations in Appendix 
[B] The analytic parameters as well as /s belong to both accretor 
masses. 

gular momentum, optimized for system age (time since the onset 
of mass transfer) f e [1 Myr, 20 Gyr], 



Md * fi-f' (A.l) 
^oi-b ^ fi ■ 
Md ^ /3 ■ 

Ma - /4-Cb 
Ma - /5-Cb 
Md ^ /6 ■ M2 

with t in yr, the mass transfer rate Md in Moyr-', the orbital pe- 
riod Porb in min and the donor mass Md in Mq. The correspond- 
ing parameters /; and for a ccretor masses M^ = 1 .4 Mq and 
10 Mq are listed in Table . 



A.l 



The exponent in the Md(f) relation 
(ei) must be roughly 1 lower than in the Md(f) relation (63), but 
not necessarily precisely 1 because both tracks have been fitted 
separately. 

For an initially 1 .4 Mq neutron star accretor, these relations 
fit the modeled helium white dwarf tracks (shown in Figs. [T[[8] 
and|9]l with an average error of 4.3%, in the tracks t - 1 Myr 
corresponds to Porb =9.0 min, Md = 2.1 ■ 10-*^ MQyr-' (lower 
than the Eddington limit; mass transfer is conservative for the 
whole range) and Md - 0.074 Mq. 

Similarly for an initially 10 Mq black hole accretor, where 
the average error is 5.3%. For a 10 Mq black hole, t - 1 Myr 
corresponds to Poib - 12.3 min, Md = 1-4 ■ 10-** MQyr-' and 
Md = 0.050 Mq. 

From the fitted parameters as well as in Fig. [T] it can be seen 
that all slopes are nearly independent of accretor mass. Also, the 
MdiPorh) relations are rather similar for both accretor masses, as 
expected. 

Appendix B: Approximate analytic solution 

By making a few mild assumptions it is possible to solve almost 
the entire evolution (assuming full feedback of angular momen- 
tum to the orbit) of a binary containing a zero-temperature white 
dwarf donor analytically, except for the early stages when the 
donor mass is high. These derivations only serve to gain physi- 
cal insight into the evolution of UCXBs - throughout the paper 
numerical solutions of more precise equations have been used. 



First we approximate the low-mass white dwarf mass-radius 
relation from Eq. (25 1 by the non-relativistic relation 



Rd/Ro = 10-2(Md/MQ) 



-1/3 



(B.l) 



which is accurate around Md - 0.01 Mq (and within 20% of 
the radius given by Eq. (Esl for 4 ■ 10-^ < Md/Mg < 1), but 



still represents the largest simplification made here. By substi- 
tuting the white dwarf radius from Eq. ( B. 1 1 into the Roche-lobe 



radius in the well-known period-mean density relation Porb = 
97r ^R^/{2GMd) we find the orbital period of a system with a 
white dwarf donor. 



Porb = AM,' 

with the constant 



A = 97:10" 



\MoK 
2G 



0.53 Momin. 



(B.2) 



(B.3) 



The orbital period is inversely proportional to donor mass (and 
proportional to donor volume). The mass transfer rate is given 
by Eq. (j6]l, in which ^d = -1/3 follows from Eq. (B.li and 
= -5/3 is the low mass ratio limit of Eq. ( [T3] l for t] - 0. 
Furthermore, we can use the gravitational wave equation and 
Kepler's third law to replace a by Porb to obtain 



Mh = -- 



48G^/3 M, 



2n 



8/3 



(B.4) 



Lastly, Md <k Ma implies Mtot ~ Ma, and we substitute Eq. (B.2 1 
for Porb to express the mass transfer rate in only the component 
masses 



Md = -CMI'^m\''^^ 
with the constant 
48 G^/^ /27rj^^^ 



C 



0.0057 M 



-'^^/3yr-' 



(B.5) 



(B.6) 



Equation (B.5i shows that the mass transfer rate depends 
strongly on donor mass. 8/3 out of the 14/3 power is due to mas- 
sive donors having a shorter semi-major axis when they fill their 
Roche lobe, which results in stronger gravitational wave radia- 
tion. The direct effect of the high donor mass is slightly smaller 
with 2 out of 14/3. Alternatively, the mass transfer rate can be 
expressed as a function of accretor mass and orbital period 



Md = -CM2/3(^) 



-14/3 



(B.7) 



By integrating the mass transfer rate given by Eq. (B.5 1, we ob- 
tain the system age f as a function of component masses. 



Ma (Im;, 



11c ' ^ 



-11/3 



(B.8) 



where integrating from 00 is justified by the negligibly short 
timescale associated with a high donor mass, as evidenced by 
the -11/3 power of the donor mass. The donor mass as a func- 
tion of time follows. 



Md 



3 

Tie 



3/11 



M-2/"f-Vl 



(B.9) 
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and by substituting Eq. (|B.2| 



'lie 



3/11 



(B.IO) 



11/3 

The evolutionary timescale is oc , even steeper than the 

gravitational wave timescale in a detached binary (oc P^^^) due 
to the effect of the donor mass. Finally, differentiating Eq. ( |B.9| l 
to time yields 



'A 



14/11 



c- 



-3/11 ^-2/11^-14/11^ 



(B.ll) 



As mentioned in Sect. 4.2.2 it shows that when comparing for 



systems with the same age, a less massive accretor corresponds 
to a higher donor mass and mass transfer rate. The mass transfer 
timescale of systems of fixed age t is given by |Md/M(j| = yf. 
This value is 2/3 times the gravitational wave timescale, and is 
independent of accretor mass in this approximation. As we saw 
in Sect. 4.2.2 however, this timescale is very slightly (~ 7% for a 
factor 7 in accretor mass) shorter for high-mass accretor systems 
when we use the more precise mass-radius relation. 

Above relations are vahd for < 0.2 Mq. Table |A.l| Hsts 
numerical values for the parameters in several equations given 
above. 



The next two expressions are only valid for q < I. The dis- 
tance between LI and the least massive binary component (2), 
valid within 0.02% (again, this is the relative error in the dis- 
tance to either component) for 7 ■ 10"^ < 1 - M2/M1 < 1, is 
given by 



c/(Ll,2) /55 
49 



+ q 



-31/42 jjj/^ ^ ^^29/72 



47 



(C.2) 



A simpler but less accurate (< 0.5%) expression, inspired by 
the form of the Eggleton Roche lobe equation ( 12 1, ia^ 



d(Ll,2) _ / ln(l+q^'^) 



(C.3) 



wh ere <q = M2 /M1 < 1. 

Warner ( 2003[ chap. 2) lists three approximations of the LI 



location, but those are valid over a much smaller range in mass 
ratio and have much larger errors than our equations. 

Appendix D: Roche lobe potential approximations 

Two simpler but less accurate versions of Eq. ([TOjl are given here. 
For a binary with Q < q = M2/M1 < 1, the potential at the Roche 
lobe surface is approximated to within 0.4% by 



B.1. Luminosity 

By assuming an accretor radius and taking the time-averaged 
bolometric luminosity of an UCXB as 



2R, 



(B.12) 



-O, 



2^1/3 



3/2 + (lOgO 

GMi/a ~ 1 - (^/5)3/2 

and to within 1 % by 

-Ol 3 



GMi/a 2 



+ 2.46q" 



(D.l) 



(D.2) 



where Md is still defined as negative, and by using Eq. (B.7i, 
we can express the following quantities in the observables L and 



and 



3/5 



(^) 



14/5 



-14/5 



(B.13) 



(B.14) 



By combining Eqs. ( B. 10 1 and ( B. 13 1 the system age follows 



(B.15) 



Appendix C: First Lagrangian point location 
approximations 

The distance between the first Lagrangian point (LI) and binary 
component 2 is given by 

d{Ll,2) ^ ^^3tanh(ln(?)/11.8)^^-0-3295-l/(21.3+0.5591n^W)J-l 

where the relative errors in c/(Ll,2) and ii(Ll, 1) are less than 
0.015% for 7 ■ 10""^ <q = M2/M1 < 1.5 ■ 10^ and less than 0.1% 
for 10"^ <q< 10^ 



" By adjusting the first term inside brackets (the 1) to 0.995, the max- 
imum error becomes 0.3%, but the approximation at ^ = 1 becomes 
worse and matches less well with its continuation at ^ > 1 . 
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